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Abstract 

^ ^. The time dependence of the survival probabihty of an opinion in a closed community has been 

^ ' investigated in accordance with social temperature by using the Kawasaki-exchange dynamics 

Q , based on previous study in Rcf. [1]. It is shown that the survival probability of opinion decays 

C/3 ' with stretched exponential law consistent with previous static model. However, the crossover 

, regime in the decay of the survival probability has been observed in this dynamic model unlike 

^ ^ \ previous model. The decay characteristics of both two regimes obey to stretched exponential. 

• Keywords: Ising Model; Politics; Random Walk; Sociophysics; Sznajd Model. 
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1 Introduction 



Binary models like Ising-type simulation have a long history. They have been applied by 
Schelling to describe the ghetto formation in the inner cities of the USA, i.e., to study phase 
separation between black and white In the sociophysics context, recently, many social 
phenomena such as election, propagation of information, predicting features of traffic, migration, 
opinion dynamics and formation in a social group have been successful modelled based on Ising 
spin systems using models and tools of statistical physics. With this respect, particularly 
successful models have been developed by Sznajd Pj, Deffuant et al.|4' and Hegselmann and 
Krause 5 . 

Among those three models, the one developed by Sznajd is the most appropriate for sim- 
ulation in networks and lattices, since it consider just the interactions between the nearest 
neighbors. Indeed, the Sznajd model has been successfully applied to model sociophysical and 
economic systems [H]- On the other hand, several modifications of the Sznajd model have been 
studied using different rules or topologies starting from different initial opinion densities [SIEIIH]- 
All these models are static (i.e. not dynamic) and they allow for consensus (one final opinion), 
polarization (two final opinion), and fragmentation (more than two final opinions), depending 
on how tolerant people are to different opinions. 

More recently the striking sociophysical model has been suggested by Aydmer ^ in order to 
explain the time evolution of resistance probability of a closed community in a one-dimensional 
Sznajd like model based on Ising spin system. It has been shown that resistance probability in 
this model decay as a stretched exponential with time. In that model spins does not move on the 
lattice sites during the simulation, so this model was so-called static. However, in a realistic case, 
spins i.e., people move in the community i.e., in the space. Social or opinion formation formed 
depend upon dynamics of the system. Because, there must be a direct connection between 
opinion dynamics and formation in a social system since the social formation is determined 
by the dynamics. Meyer-Ortmanns studied recent work in which the condition for ghetto 
formation in a population with natives and immigrants by using Kawasaki-exchange dynamics 
in a two dimensional Ising model. She showed that ghetto formation can be avoided with a 
temperature increasing with time. Similarly, Schulze have also generalized Meyer-Ortmanns 
work to up to seven different ethnic groups to explain ghetto formation in a multi-cultural 
societies in a Potts- like model [TU) . 

In this study, we have developed a dynamic version of the Aydmer (Tl model by combining 
the Aydmer and Meyer-Ortmanns [Hj models based on one-dimensional Ising model. 

2 Kinetic Model and Simulation 

In one-dimensional static model m, each site carriers a spin which is cither spin up (-1-1) or 
spin down (-1) randomly. Spin up (-1-1) represent the host people and spin down (-1) represent 
the soldier. The host people always against occupation, and, on the other hand, soldier always 
willing to continue occupation, who always have the opinion opposite of that of the host people. 
Furthermore, the comnmnity member i.e., spins doesn't also move on the lattice during the 
process. 

In this model, initially, it was assumed that there was a over all consensus among member of 
the community against occupation even if some exceptions exist. One expects that host people 
obey to this consensus at least initially. In this sense, community behaves as polarized at zero 
social temperature \1'6\ against occupation just like Ising ferromagnet at zero temperature. 

It was conjectured that host people are infiuenced by soldiers even though they against 
occupation owing to they are exposed to intensive biased information or propagation. Soldiers 
affect the host people and force to change their opinion about occupation. Effected people may 
change their own opinions depending on resistance probability of the nearest neighbors about 
occupation. Moreover, effected host people affect neighbors. Such a mechanism depolarize the 
polarization (resistance probability) of all host people. Hence social polarization destroy. 
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However, soldiers, unlike host people, have not been influenced by the host people. Their 
opinion about justifying the occupation does not change during the occupation process, since 
they may be stubborn, stable or professional etc., who behaves like persistent spins in Ising 
spin system. It is means that the probability of the against occupation of a soldier is always 
zero. 

If we summarize, we can say that none spins does flip fully in the system. Spin up always 
remains spin up, and spin down always remains spin down. In this respect, the probability of 
against occupation of host people can be interpreted as a survival probability of opinion of host 
people about occupation under above considerations. In this sense, the survival probability 
Wi of opinion of host people indicate equal to 1 at least initially and, on the other hand, the 
probability of against occupation of soldier equal to zero, which means that soldier behaves as 
a trap point lattice which depolarize the survival probability of opinion of host people. 

Of course, one may suggest that there are many different number of opinions in society, 
however, it is possible to flnd that a society being formed two-state opinion in a real case. 
Therefore this model is a good example for two-state opinion model as well Galam contrarian 
model even though it seems that it is very simple. Furthermore, in real social systems, 
people move on the space, i.e., lattice. Therefore, in this study, we assumed that people 
i.e., spins randomly move on the lattice through the Kawasaki-exchange dynamics contrary to 
previous model. 

The survival probability Wi for a people at site i at the next time t -f 1 is determined with 
the survival probability of nearest-neighbors with previous time t as 

W,it+l)^^[W,+i{t) + W,-i{t)]. (1) 

We note that the survival probability for all site are calculated as synchronously. 

Randomly motion of the spins i.e., people on the lattice through the Kawasaki-exchange 
dynamics. Firstly, a spin pair is chosen randomly and then it is decided whether spin pair 
exchange with each other or not. In this approach, the nearest-neighbor spins are exchanged 
under heat-bath dynamics, i.e., with probability p ~ exp (—AE/ksT)), where AE is the energy 
change under the spin exchange, ks is the Boltzmann constant, and T is the temperature i.e., 
social temperature or tolerance. Hence, to obtain probability p we need to calculate Ei and E2 
which correspond to energy of the spin pair at first position and after exchange with position 
of spins, respectively. Energy Ei and E2 can be calculated in terms of the survival probability 
instead of spin value as 

Eiit) = aW^it)+bW,+iit) (2a) 



E2{t) = aW,+i{t) + bW^{t) (2b) 

where 

a=[W,-iit) + W,+iit)] 

and 

b=[W.dt) + W,+2{t)] 

Energy difference is written as AE ~ E2 — Ei from Eq. (Pa|l and (|2b|l . 

In addition, the total survival probability of opinion of host people at the any time t can be 
obtained over each person for any r configuration as 

Prit) = —J2W,it) (3) 

where toq is the initial number of host people. On the other hand, the averaged survival 
probability at the any time t can be obtained from Eq. |(2Jl over the independent configuration 
as ^ 



r— 1 



where R is the number of different configurations. 
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Fig. 1: Data for p = 0.1 were plotted In < P{t) > versus t in Fig. Q^a), log < P{t) > versus 
logt in Fig. mb), < P{t) > versus Int in Fig. Q^c), and log(— In < P{t) >) versus \ogt in 
Fig. md). Solid-dot lines indicate data and solid lines represent fitting curves in all figures. 

3 Results and Discussion 

We have adopted the Monte Carlo simulation technique to the onc-dimcnsional sociophysical 
model using the lattice size L — 1000 with periodic boundary condition, and independent 
configuration R = 1000 for the averaged results. The simple algorithm for the simulation is 
as follows: i) at the t = 0, Eq. Q is initially calculated, ii) for t > a spin pair is randomly 
chosen, and then it is decided whether the spin pair exchange or not with the probability 
p ~ exp (— AE'/fcsT)), this step is repeated L times, iii) after ii-steps are completed, Eq. Q is 
recalculated again, and to continue this procedure goes to step ii. 

The simulation results are as follow: We have firstly plotted simulation data versus time in 
Fig. n in a several manner. It is explicitly seen from Figs, ^a)-(c) that there are no power, 
exponential and logarithmic law dependence in our simulation data, respectively. However, as 
seen Fig.^d), data well fit to the stretched exponential function as 

< P{t) >^ e-^*" (5) 

where A is the relaxation constant, and (3 is the decay exponent of the survival probability. 
This result indicate that the time evaluation of survival probability of the opinion of the host 
people in a closed community has stretched exponential character i.e., Kohlraush- William- Watts 
(KWW) decay law [111112 • 

It should be tested whether Fig.^d) satisfies to stretched exponential or not |^. Because, 
as noted by Stauffer, the Fig.QId) would work as stretched exponential, if pre-factor of Eq. ^ 
is equal to 1. However, if pre-factor is less than 1, it may give the impression of stretched 
exponential form, even for (3 — 1. Therefore, it can be plotted — In < P{t) > versus suitable 
powers of t, like i^/^, t^/"^, etc., and find out the best straight line among the powers of t for long 
times. Hence, — In < P{t) > was plotted versus powers of t for p = 0.1 then the best straight 
fitting hue for long times was obtained for /3 = 0.8 for T = 0.0001, (i = 0.7 for T = 0.01, and 
(3 — 0.6 for T = 0.1 as seen in Fig. |nia)-(c) respectively. These results confirm to this method 
used to find out stretched exponential exponents in Fig. Q^d), and also all figures in Fig. [3 as 
mentioned . Also, this test indicates that prcfactor in Eq. (jSJ does not effect results presented 
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Fig. 2: The time dependence of the survival probability of the opinion of host people de- 
cays KWW i.e., stretched exponential with time for different social temperature T. The time 
crossover appears in the time evolution of survival probability of the opinion at low social 
temperatures. The crossover becomes more clear when social temperature decreases. 




Fig. 3: (a) - In < P{t) > versus t°-^ for T = 0.0001, (b) - In < P{t) > versus f-'^ for T = 0.01, 
and (c) — In < P{t) > versus for T — 0.1. All figures are plotted for fixed p value i.e. 
p = 0.1, and solid lines represent fitting curves. 



5 




Fig. 4: (a) The soldier density dependence of the exponent f3 under and below the crossover 
time tc for different social temperatures, (b) Change of the decay exponent (3 depend on social 
temperature T under and below of the crossover time tc for a fixed soldier density p = 0.1. 



in this paper. 

It is concluded that results for high temperatures also consistent with static model Jj. But, 
unlike the static model, time crossover has been observed in dynamic model at low tempera- 
tures. In order to investigate the transition we have plotted survival probability versus time for 
different social temperature T in Fig. [3 It is clearly seen that the time crossover occurs depend 
on social temperature. When social temperature decreases, the crossover become more clear. 
Such a behavior was not observed in a static model. We can bridge the short time regime and 
the long time regime by a scaling function / {t/tc) 

(P(i)) = e-(*/-)V(iAc) (6) 

where tc indicates the time crossover. For our simulation data, the scaling relation (6) can be 
written for very long and very short time intervals as 

r p-it/rf^ if t«t 



On the other hand, in order see how the decay exponent [3 depend on soldier density p, and 
social temperature T, we have plotted (5 versus soldier density p in Fig. ^a) for t < tc and 
t > tc in account to taken different social temperatures, and social temperature T in Fig. ^b) 
for a fixed value of density p, respectively. 

As seen from Fig. EJa) that /3i and (32 are linearly depend on soldier density both of two 
regimes at low social temperature. On the other hand, the decay exponent has two different 
character for t < tc and t > tc depend on social temperature T in Fig.^fb), the decay exponent 
Pi decreases with increasing temperature T for t < whereas /32 increases with increasing 
temperature T for t > tc 'At low temperatures. However, for relatively high temperatures we 
roughly say that /3i approach to /32 for both two regimes obey to Eq. |(7|). 

Finally, to understand the social temperature and soldier density dependence of the time 
crossover t^ we have plotted tc versus social temperature T in Fig. Ela) for a fixed soldier 
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Fig. 5: (a) The social temperature dependence of the crossover time tc for fixed soldier density 
p = 0.1, (b) The soldier density dependence of the crossover time tc for fixed social temperature 
T = 0.0001. 

density p, and versus soldier density in Fig. El^b) for fixed social temperature T, respectively. 
It seems from Fig. I^a) that the crossover transition tc quite rapidly decrease with increasing 
T, on the other hand, it seems from Fig. ]^h) that it slowly decrease with increasing soldier 
density p. We note that as seen inserted figure in Fig. IHb) the crossover transition tc depends 
on soldier density with power law for fixed social temperature. 

4 Conclusions 

We suggest that the stretched exponential behavior of decay must be originated from model 
system. The persistent spins i.e., the soldiers doesn't flip during simulation, therefore they 
behave as a trap in the system. Hence they play a role diminishing the survival probability 
of the neighbor spins in the system. Consequently, decay characteristic of the system can be 
explain due to the trapping states. Another say, this characteristic behavior doesn't depend on 
either diffusion dynamics of spins or interaction rules between spins. 

Another unexpected behavior is the time crossover in P contrast to previous model [I] . We 
supposed that this amazing result originated from opinion dynamics depend on social temper- 
ature. Model allows to the opinion formation with time. Indeed, there is a direct connection 
between opinion dynamics and formation in a social system since the social formation is de- 
termined by the dynamics as depend on the social temperature. For example, in a real spin 
system, decreasing temperature phase separation may occur in the system. In the sociophysical 
sense, it means that people who have different opinion are separated each other with decreasing 
social tolerance, and therefore the ghetto formation or polarization may occur in the system. 

It is expected that interactions between soldier and host people is maximum when soldiers 
are randomly distributed in the community. As social temperature, i.e., tolerance is decreased, 
however, phase separation occur with time, so this leads to decreasing of the interactions. 

In our opinion, the ghetto formation in the system doesn't leads crossover transition in 
time because of the ghetto formation is randomly distributed relatively. On the other hand, 
the time average of survival probability over different configuration effect of ghetto formation 
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may probably destroy. So we don't hope that ghetto formation is not responsible crossover 
transition. However, polarization must be occurred at low temperature leads to meaningful 
phase separation in the system. Such a polarization may leads to crossover transition in time. 

Stretched exponential behavior indicates mathematically that decay for the relatively short 
times is fast, but for relatively long times it is slower. One can observe that this mathematical 
behavior corresponds to occupation processes in the real world. In generally, a military occu- 
pation is realized after a hot war. The community does not react to occupation since it occurs 
as a result of defeat. People are affected easily by propaganda or other similar ways. Therefore, 
it is not surprised that resistance probability decrease rapidly at relatively short times. On 
the other hand, spontaneous reaction may begin against occupation in the community after 
the shock. Hence, community begins by regaining consciousness and more organized resistance 
may display difficulties for occupants. For long times, the resistance probability decreases more 
slowly. This means that resistance against occupation extends to long times in practice. At 
this point, the number of soldiers is also important, because the density of soldiers determines 
the speed of decaying. 

The different regimes have been observed in the decay of the survival probability. These 
regimes clearly appear particularly at low temperatures. In the case of the social temperature is 
very low. Pi is bigger than /32 which indicates the decay of the survival probability for relatively 
short time is slower than for relatively long time. This can be interpreted that the resistance of 
host people against occupation may be broken spontaneously if soldier can wait enough time. 

Of course, the mechanism considered in this work can be regarded as simple, but, it would be 
useful to understand the time evolution of the resistance probability of the community against 
to occupation in the one-dimensional model under some considerations. We remember that 
simple social rules lead to complicated social results, hence we believe that the obtained results 
and model can be applied the real social phenomena in the societies to understand the basis of 
them. 
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